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ABSTRACT 

We report results of a systematic study of the broad band (2-2000 keV) 
time resolved prompt emission spectra of a sample of gamma-ray bursts (GRBs) 
detected with both Wide Field Cameras on board the BeppoSAX satellite and 
the BATSE experiment on board CGRO. In this first paper, we study the time- 
resolved dependence of the intrinsic peak energy E p>i of the EF(E) spectrum 
on the corresponding isotropic bolometric luminosity L iso . The E Pti -L- lso relation 
or the equivalent relation between E pi and the bolometric released energy E iso , 
derived using the time averaged spectra of long GRBs with known redshift, is 
well established, but its physical origin is still a subject of discussion. In addition, 
some authors maintain that these relations are the result of instrumental selection 
effects. We find that not only a relation between the measured peak energy E p 
and the corresponding energy flux, but also a strong E Pii versus L iso correlation 
are found within each burst and merging together the time resolved data points 
from different GRBs. We do not expect significant instrumental selection effects 
that can affect the obtained results, apart from the fact that the GRBs in our 
sample are sufficiently bright to perform a time-resolved spectroscopy and that 
they have known redshift. If the fundamental physical process that gives rise to 
the GRB phenomenon does not depend on its brightness, we conclude that the 
found E Pt i versus Li SO correlation within each GRB is intrinsic to the emission 
process, and that the correlations discovered by Amati et al. and Yonetoku et 
al. are likely not the result of selection effects. We also discuss the properties of 
the correlations found. 

Subject headings: gamma ray bursts: general — gamma rays: stars 
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Introduction 



In spite of the major advances in knowledge of gamma-ray bursts (GRBs) a fterglow 
prop erties, mainly made with Swift, the GRB phenomenon is still poorly understood (ILyutikov 
20091 ). It is recognized that the study of the prompt emission is of crucial importance, as 
it is more directly connected with the original explosion. One of the issues still open is the 
radiation emission mechanism(s) at work. Most of the GRB properties derived thus far come 
from the time-averaged spectra. The function that has been found to be tter describe them 



from 15 keV up to 10 MeV is a smoothly broken power-law proposed by iBand et al.l ( 11993 
Band function, BF). On the basis of the spectral data obtained with the Burst and Tran- 
sient Source Experiment (BATSE), aboard the Compton Gamma Ray Observatory sat ellite 
(CGRO) and, for example, with the BeppoSAX GRBM data (e.g.. iGuidorzi et al.ll201lh . for 
long GRBs (>2 s), the mean value of the low-energy photon in dex a of the BF is about 
—1, while that of the high-energy photon index /3 is about —2.3 (IKaneko et al.ll2006l ). As a 
consequence of this result, the received power per unit logarithmic energy interval EF(E) 
shows a peak value that, in the BATSE era, seemed to show a sharp Gaussian distribution 
around 200 keV. With the discovery of the X-ray flashes with BeppoSAX , later also found 
with HETE-2, Swift, and, now, wit h the Fermi Gamma-ray Burst Monitor, th is distribution 
results in being much flatter (e.g., iKippen et al.ll2003l ; ISakamoto et al.ll2005l ). In the cases 
in which j3 cannot be constrained, a power-law model with a high energy exponential cutoff 
(CPL) gives a good fit to the data, and, in some cases, even a simple power-law can descr ibe 
the GRB time averaged spectra up to several MeV photon energies ( IKaneko et al.l 120061 ). 



Several radiative emission models have been worked out for the interpretation of the 
GRB spectra. Given their non-thermal shape, the first mo del proposed was synch r otron 
emission by non t hermal electrons in strong magnetic fields (IRees fc Meszarod Il994j ; iKatz 
1994 ; lTavanilll996l ). Indeed, the time- averag ed spectra of many GRBs are c onsistent with an 
optically thin synchrotron shock model (e.g.. lTavanilll996l ; lAmati et al.ll200ll ). However, there 
are a significant number of GRBs for which this model does not work. Indeed, while for an 
optically thin synchrotron spectrum, the expected power-law index of the EF(E) spectrum 
below the peak energy E p cannot be steeper than 4/3 (ideal case of an insta ntaneous spectrum 



Preece et al 



1998 



in wh ich electron cooling is not taken into account), in many cases (e.g.,_ 
2000[) the measured spectra, even those time resolved (ICrider et al.l Il997t iFrontera et al.l 
2000), are inconsistent with these expectations. 
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of these models interprets some of the prompt emission features, but fails to interpret others. 



One of the GRB spectral properties that the emission models should interpret is the 
correlation between the intrinsic (redshift c orrected) peak ene rgy E Pii of the EF(E) function 
and either t he GRB released ener gy E- lso (jAmati et al.l 120021 ) or the peak bolometric lumi- 
nosity -Lp i SO (jYonetoku et al.ll2004j ). Both correlations (the Amati relation and the Yonetoku 
relation) have been derived from the time-integrated spectra assuming isotropic emission. 
The Yonetoku relation followed the Amati's, from which the so-called Ghirlanda relation 



(IGhirlanda et al.ll2004j ) was also derived by replacing the released energy E iso with that (£" 7 ) 
corrected for the beaming factor (E 1 = (1 — cos^)^^). The latter is model dependent, 
being derived assuming a jet like structure of the fireball, a constant efficiency of the fire- 
ball in converting kinetic energy in the ejecta into gamma rays, a mass density distribution 
of the circumburst medium, and, mainly, that the break time observed in the late after- 
glow light curve occurs when reciprocal of the bulk Lorentz factor of the jet, 1/r, becomes 
of the order of the jet opening angle #j et . The latter assumption re quires that the break 



time is achromatic, p roperty not observed in many Swift GRBs (e.g.. iCampana et al.ll2007l ; 
Melandri et al.l 120081 ). In any case, these relations are all equivalent as far as the physics: 
a relation between the photon energy in which most of the energy is released and the elec- 
tromagnetic radiation emitted by a GRB, which can be expressed equivalently in terms of 
released energy or peak (or average) GRB luminosity. 



The Amati relation (jAmati et al.l 120081 . E, 



0.54 ± 0.03, where E Pti is measured in keV and E is 



KE™ Q , with K = 98 ± 7 and m = 
is given in units of 10 52 erg), is satis- 
fied, within an extra-Poissonian scatter of log E p ^ normally distributed around the best-fit 
power-law with o ~ 0.2 dex, by all long GRBs (more than 100) with known redshift z dis- 
covered thus far, except the nearest and leas t energetic GRB ( z = 0.0085) ever observed 
(GRB 980425) and, maybe, GRB 031203 (e.g. kmati et alil2007h . but not by short GRBs. 
Also the Yonetoku re l ation is satisfied for GRB s with known redshift, except GRB 980425 



(e.g., 

L, 



p,iso and E\ so 



(IGhirlanda et al 



2005bJ). 



Ghirlanda et al.l (2005b); INava et al.l ( 120121 )). given the tight correlation found between 



In spite of this, the Amati relation has been questioned by various authors ( 


3and & Preece 


2005 




Butler et al. 


2007 




2009 




Shahmoradi & Nemirofi 


2009 




Collazzi et al. 


2011 


), main- 



taining that it is likely the result of selection effects , even if these effects, when investi- 
gated by other authors (IGhirlanda et all l2005al . I2OO8I : iNava et al.l boosl : kmati et all 12009 



Krimm et al.ll2009l ; INava et al.l 1201 ll ) are found to be marginal. 



One of the major criticisms of the A mati relation is that its norm alization depends on 
instrument sensitivity (IButler et al.l 120071 ) . However, this conclusion by iButler et al.l (120071 ) 



for Swift GRBs was not based on measured spectral peak energies E P)i , but on E Pt i values 
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inferred, under som e assumptions, from a Bayesian method. An investigation, performed by 
Amati et al.l (120091 ) by deriving the E p>i dependence on E iso for different sets of GRBs, each 



obtained with a different instrum ent, and by usin g the m easured values of E Pti reported in 
th e Swift BAT official catalog by ISakamoto et al.l ( 120081 ). has not confirmed the inferences 
bv lButler etaD (j2007h . 

Other works that question the reliability of the E P)i - E iso make use of the (E p , Fluence) 
observer plane, positioning in this plane the data points obtained from the spectral analy- 
sis of GRBs detec t ed wi th different satellite instrum ents, using the method first proposed 



by iNakar &: Piranl (120051 ). While some authors (e.g. iBand &: Preecd l2005t iGoldstein et al. 



201qT T find that most of the l ong BATSE GRBs do not satisfy the relation, other authors 



(e.g., iGhirlanda et al.l l2005at iNava et al.l 120081 1201 ll ) find that most of them do. This dis- 
crepancy is due to various reasons, such as the condition assumed to satisfy or not satisfy 
the Amati relation and the systematic errors in the determination of fluence and E p (e.g. 



Collazzi et al.l 1201 ll ). However some results are well established and are shared by most au- 
thors who have performed the E P j - Ei SO test. GRBs with known redshift are not a special 
sub-pop ulation, they are even l y distributed along the weaker correlation between E p and 
fluence (IGhirlanda et al.ll2008l : fcollazzi et aJboiA. No evidence of evolution of th e E p ^ - 
Ei SO correlation with the redshift is found ( IGhirlanda et al.ll2008l ; INava et al.l 120 111 ). In the 
observer frame, after all the results obtained with Swift and Fermi, no bursts with large 
fluence and low/intermediate E p have been found. They do not exist or should be very rare: 
Nothing prevents their detection. Instead, bursts with intermediate/high E p and small flu- 
ence could be affected by instrumental selection effects, such as the minimum flux to trigger 
a GRB, the minimum fluence to fit its spectrum and constrain its E p , and truncation effects 
related to the instrument passband (ILloyd et al.ll2000l ) . The case of Swift BAT is an example 
of instrument affected by truncation biases: In spite of its higher fluence sensitivity than 
BATSE, E p can be accurately determined only if it has a value within or very close to its 
energy passband (15-150 keV). However, we find so me results abou t the v alidity of the E Pti 
- E[ SO relation untenable, such as the conclusion by I Collazzi et al.l ( 120121 ). who find that a 
significant fraction of the same GRBs with known redshift which have been used to derive 
the Amati relation do not satisfy their test. Clearly, their test condition (their so-called 
"Amati limit") in the (Fluence, E p ) plane is too restrictive. The same authors state that, 
while the Amati relation is the result of selection effects, the Ghirlanda relation is valid. 
Also, this statement is problematic, given that this relation is based on the same E p and 
Fluence measurements, being the only correction performed with, as discussed above, the 
replacement of E iso with E T In addition, the "Ghirlanda limit" is derived using rather loose 
assumptions about the beaming angle. 



The opposite conclusion was recently reached by 



Yonetoku et al 



(120101 ) via analyzing in 
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detail all possible data truncation and detector sensitivity effects, and by lNava et al.l ( 120121 ). 
who analyzed a complete sample, for redshift determination, of bright Swift GRBs (1 s peak 
photon flux P > 2.6 photons s _1 cm -2 in the 15-150 keV BAT band). 

From this long standing debate, it is apparent the need to explore other approaches 
for testing the origin of the spectrum-energy correlations, i.e., whether they are the re- 
sult of instrumental selection effects or are related to the fundamental physics of the GRB 
phenomenon. Given the significant evolution of the GRB spectra, studies of time-resolved 
spectra are crucial not only to test the E Pti versus L iso relation, but also to delve deeper into 
the issue of the radiative mechanisms at work during the prompt emission. 

Motivated by both these needs, we performed a systematic study of the broad band (2- 
2000 keV) time-resolved prompt emission spectra of a sample of GRBs detected with both 
Wide Field Cameras (WFCs) aboard the BeppoSAX satellite and the BATSE experiment 
aboard the CGRO. In this paper we will concentrate on the test of the E P j versus L iso relation. 
The WFCs were among the few instruments that detected GRB prompt emission down to 
2 keV. Thus we can obtain time-resolved spectra in an energy band still not well explored: 
the 2-2000 keV band. A paper devoted to testing physical emission models of GRBs using 
the same time-resolved spectra will be the subject of a forthcoming paper (Frontera et al. 
2012, in preparation). 



2. The GRB sample and spectral analysis 



There were nine GRBs simultaneously detected with WFCs and BATSE: 970111, 971206, 
971214, 980329, 980519, 990123, 990510, 990907, 991030. We performed an analysis of the 
prompt emission spectra of all of them. However, only four of them (970111, 980329, 990123, 
and 990510) were sufficiently bright to allow a fine time resolved spectroscopy (see Tabled]). 

The instrumentation that dete cted these burst s is wi dely described in the lit erature. For 
the B ATSE experiment, see, e.g., iFishman et al.l (119941 ). while, for WFC, see iJager et al. 



(119971 ). The BATSE spectra were taken from the Large Area Detectors (LADs), whose 
typical passband is 25-2000 keV. The LADs provide various types of data products. The 
data used for this analysis are the high energy resolution burst data, that provide 128 energy 
channels with a minimum integration time of 64 ms. For details of the BATSE spectral data 
products and detector response matrix see iKaneko et al.l (120061 . and references therein). 



The WFCs consisted of two coded aperture cameras, each with a field of view of 40° x 40° 
(full width at zero response) and an angular resolution of 5 arcmin. They operated in normal 
mode with 31 energy channels in 2-28 keV and 0.5 ms time resolution. 
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The background-subtracted light curves of the four strongest GRBs in our sample, 
detected with both BATSE and WFCs, are shown in Figure HJ For the BATSE data, 
the background level was estimated using the count rates immediately before and after the 
GRBs. Given that the background is variable during the GRB, it was estimated by means 
of a parabolic interpolation, channel by channel, between the background measured before 
the event and that measured after the event. For WFC spectra, the background level was 
estimated using an equivalent section of the detector area not illuminated by the burst or by 
other known X-ray sources. We also checked the consistency of this background level with 
that obtained by using the data before and after the burst. 

We subdivided the time profile of each GRB into a number of time slices (see Figured]), 
taking into account the GRB profile as observed with 1 s time resolution (visible pulses, their 
rise, peak, and decay) and the count statistics. We performed further spectral analysis in all 
of the time slices in which it was possible to constrain the peak energy E p . With reference to 
Figure[U in the case of GRB 970111 we excluded the first two (nos. 1 and 2) time slices for 
their low count statistics. For the same reason, we excluded from the analysis the first (no. 
1) and the last (no. 8) of GRB 980329, the first (no. 1) of GRB 990123, and intervals 3, 4, 
5, and 6 of GRB 990510. For this GRB, we also excluded intervals 13, 14 and 15, given that 
only an upper limit could be obtained for E p . In the case of GRB 990123, we also excluded 
from the analysis the time intervals from 21 to 26, given that this part of the event was 
observed by WFC through the Earth's atmosphere. The number of time intervals in which 
we subdivided the time profile of each event, the number of selected time slices in which it 
was possible to estimate and constrain E p , the GRB fluence, and its redshift when known, 
are given in Table [H 

For each WFC+ BATSE time resolved spectrum we used as an input model the BF: 




if 



(a -(3) 



E > E 



and 




if 



OS -at) 



E < E 
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where a and /3 are the power law low-energy (below Eo) and high-energy (above Eq) 
photon indices, respectively, and A is the normalization parameter. 

In the fit, along with A, a, and (3, we adopted as a free parameter, instead of E , the 
photon peak energy E p = E (2+a). In addition, a normalization factor between BATSE and 
WFC data was included in the fit and left to vary in the range 0.8-1.2, to take into account 
a possible intercalibration error. Actually, we found that, for all analyzed GRB spectra, 
this parameter was consistent with 1. The systematic error used by the BATSE team to 
take into account the uncertainty in th e background subtrac tion and the uncertainty in the 
instrument response function (see, e.g.. iKaneko et al.l (120061 ) for details) was injected in the 
fit. 

The input model was assumed to be photoelectrically absorbed ( WABS model in XSPEC). 
Given that the absorption colu mn density Nw could not be constrained, a Galactic absorp- 
tion along the GRB direction (IDickey fc Lockmanlll990l ) was assume d. To deconv olve the 
count rate spectra we adopted the xspec (v. 12.5) software package (jArnaudlll996l ). If not 
explicitly stated, the quoted uncertainties are single parameter errors at the 90% confidence 
level. 



3. Results 



The fit results of the BF to the joint WFC + BATSE time-resolved spectra of GRBs 
970111, 980329, 990123 and 990510 are reported in Tabled The time behavior of the best-fit 
parameters to the tested model is shown in Figure [2J In the bottom panels of this figure, 
the time behavior of the null hypothesis probability (NHP) is also shown. 

From Figure [2] it is apparent that E p , for each event, mimics the time behavior of the 
2-2000 keV flux. This result is better shown in Figure |3l where the time-resolved peak 
energy derived from the best fit of BF to the time-resolved spectra is plotted as a function 
of the corresponding flux. As can be seen from this figure, a positive correlation is found for 
each of the four GRBs, consistent with a power-law dependence of the peak energy on flux, 
even if the statistical significance changes from one event to another, also depending on the 
number of available data points (see Table [3]). The power-law parameters were derived in 
two ways: using the simple least-squares method, and the maximum likelihood method in 
the case that the correlated data (xi,yi) could be described by a linear function Y = mX + q 
with the addition of an extrinsic (non Poissonian) variance a^ xt among the free parameters 



( ID'Agostini 



(see, e.g., 



2005). The latt er method has already been adopted for various applications 



Amati et al.l (120081 ) and references therein). 
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Given that the peak energy E p is related to both Eq and a (E p = Eo(2 + a)), in order to 
establish the a contribution to the strong correlation between E p and flux, we investigated 
the behavior of a with E p . The result is shown in Figure HI As can be seen, even if, in two 
cases (GRB 970111 and GRB 980329), we see a rapid change of a with E p , no significant 
correlation is found in the cases of the other GRBs, especially in the case of GRB 990123, 
where we have the highest number of time-resolved spectra. 

For GRBs with known redshift (GRB 990123 and GRB 990510), using the concordance 
cosmology (Q\ = 0.73, Q m = 0.27, H = 70 km s _1 Mpc -1 ), we derived the intrinsic time 
resolved peak energy E p ^ as a function of the corresponding 2-2000 keV isotropic luminosity 
L iso . The result is shown in Figure [5j It does not change if we integrate L iso up to 10 MeV 
using the best fit parameters of the BF. As can be seen from this figure, and quantified by 
the best fit correlation results (see Table @J, E p i is related to L iso through a power-law, with 
a very low probability that the result is due to chance, especially in the case of GRB 990123. 
For each of these GRBs, we also derived Li SO separately in the energy band below E Pt i and 
above E p ^. The result is that in both ranges the correlation between E p ^ and Li SO still holds. 

In order to see the cumulative result, we have merged all the available data on GRBs 
990123 and 990510, to derive the mean dependence of E p ^ on Li so . The result is shown in 
Figure [61 A power-law dependence of E Pii on L iso is still found (see Table @J, with a very low 
null hypothesis probability (1.57 x 10~ 13 ), much lower than that found within each event, 
confirming the reality of the E p -L lso correlation in the single events. 



Discussion and conclusions 



By joining together the WFC and BATSE spectral data of the four strongest (> 15 x 
10 -6 erg cm -2 ) GRBs simultaneously observed with both instruments, it was possible to 
perform a fine time-resolved spectral analysis in the broad energy band 2 keV to 2 MeV, a 
passband still scarcely explored SIS db whole. 

We do not expect significant systematic errors from this joint analysis. The response 
functions of both instruments are well known. Indeed, in the fits the cross-calibration factor 
was found to be always 1, in spite of being left to vary between 0.8 and 1.2 (see Sec- 
tion [2]). This is not the first time that a joint WFC/ BATSE spectral analysis has been 
perfor med. Results of s imilar analyses have bee n reported in the past by the BATSE 
team fariggs et alihoool : iKippen et al.l fiooil . 120041 ). In addition, the BATSE -deconvolved 
spectra of bright G RBs were cross-checked with those obtained with BeppoSAX GRBM 
(IFrontera et al.l 120091 ) : these in turn were cross-calibrated with WFC, with many published 
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results (e.g. iFrontera et al.lll998l . l2000f ). 



For each of the strongest GRBs, we obtained a significant number of time-resolved 
spectra with constrained E p : 8 spectra for GRB 970111, 6 spectra for GRB 980329, 19 spectra 
for GRB 990123, and 7 spectra for GRB 990510, with a total number of 40 analyzed spectra. 
With these spectra, we investigated the dependence of the time-resolved peak energy E p on 
the corresponding 2-2000 keV flux using the empirical BF as an input model. 

We find a significant power-law correlation between the derived peak energy E p and the 
flux within each GRB in our sample. The power-law best-fit parameters, evaluated with two 
different methods (the least-squares method, and the likelihood method with the addition of 
an external variance as a free parameter) give (see Table [3]) similar index values in the cases 
of GRBs 970111 (0.68 ±0.06 versus 0.65^;^) and 980329 (0.16 ±0.04 versus 0.15±g;gf), and 
different values in the case of the GRBs 990123 (0.53 ± 0.05 versus 0.46±g;g|) and 990510 
(0.81 ± 0.15 versus 0. 56^23), even if, in the latter case, these values are statistically almost 
consistent with each other. However, the likelihood method has the advantage of giving us 
information about the non-Poissonian spread of the data points around the best-fit curve. 
It shows that such a spread, even if small, affects the correlation, with the highest value 
(o"cxt = 0.09 dex) in the case of GRB 990510 and the minimum one (<T ex t = O.OO^qq dex) 
in the case of GRB 980329. Actually, in this case, also due to the low statistics of the data 
points, the correlation significance is low and, as shown by the reduced y 2 value (see Table [3]), 
is not sensitive to an external spread. 

No clear correlation is found between the low-energy photon index a and flux and 
between a and E p , apart from in one case (GRB 980329). This means that the most relevant 
parameter that gives rise to the E p -Hux correlation is the break parameter Eq in the BF. 

An equivalent important result is that the power-law correlation found between E p and 
flux is confirmed when we correlate, for two GRBs with known redshift, the intrinsic peak 
energy E p ^ with the corresponding isotropic luminosity. Notice that the power-law index of 
the correlation changes from one GRB to the next when we use the least-squares method 
(0.53 ± 0.05 for GRB 990123, and 0.81 ± 0.15 for GRB 990510), but does not when the 
likelihood method is used (0.46 ± 0.09 for GRB 990123, 0.56+°;| for GRB 990510). This 
fact clearly means that the power-law index is affected by the external spread. If this spread 
is due to an unknown physical parameter, this result is an important hint for physical models 
of the prompt emission process: The intrinsic E Pti derived from the assumed emission model 
should be related to at least two physical parameters of the model. 



When we join together all the available data on GRBs with known redshift we find that 



the Epi-Li SO correlation becomes even more robust, with a probability of 1.57 x 10 13 that 
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the correlation averaged over all data is due to chance, and a power-law index that is almost 
independent of the used best fit method (0.66 ± 0.03 in the case of the least-square method 
versus 0.63^oq7 in the case of the likelihood method). It is also interesting to note that the 
non-Poissonian spread found for the GRB averaged corr elation (cr ext = Q1Q 6-005) is three 
times lower than the spread found in the Amati relation ( lAmati et al.l 120091 ). which is well 
known to be based on time-averaged spectra. This is a strong hint that part of the spread of 
the Amati relation is related with the fact that the E p values determined from time-averaged 
spectra are biased because of the spectral evolution of the GRB prompt emission. 



If we compare our GRB-averaged correlation result with that obtained by lGhirlanda et al 



(120101 ) from the time-resolved spectra of Fermi GRBs, we find that our results are consistent 
with those within a 2a belt, even if the slope obtained by these authors (0.36 ±0.05) is lower 
than that found by us (0. 63^07; see Table 0] and Figure [6]), likely due t o the sample var iance. 



Indeed, our found slope is similar to that (0.621 ± 0.003) reported by iLu et al.l (120121 ). who 
performed the time-resolved spectral analysis of a sample of 15 Fermi GRBs with known 
redshift. 

Within a 2a spread, o ur results are also consi stent with the time-averaged E P)i versus 
Lj SO correlation, obtained by lGhirlanda et al.l (120101 ) using 95 pre- Fermi plus 10 Fermi GRBs 



(see the right panel of Figure [7]) . 

We cannot exclude that selection effects can influence our results, as with the fact that 
the analyzed GRBs are detected by both instruments, that they are sufficiently bright to 
allow time resolved analysis, and that they have known redshift. However, it seems difficult 
that this unavoidable selection can introduce a correlation between E p and L; so within single 
GRBs. In addition, a similar correla tion has been found by the other mentioned authors 
( jGhirlanda et al.ll2010t ILu et al.ll2012l ) with other GRBs and with other instruments. 



To conclude our results strongly s upport, at least in the range of lumi nosities explored 



with our data, the reality of the Amati ( lAmati et al.ll2002l ) and the Yonetoku (jYonetoku et al. 
2004 ) relations, both derived using time averaged spectra. 



Also, our results give strong constraints on the physical models. In a forthocoming paper 
(Frontera et al. 2012, in preparation), with the same data we are testing different physical 
models, among them the recently developed grb comp model, which is devoted to the spec- 
tral formation of a GRB ( jTitarchuk et al.l |2 0121 ). In this model a physical interpretation of 
the Amati relation is also given. 



We are grateful to Pawan Kumar for valuable discussions and suggestions. The Bep- 
poSAX satellite was a joint effort of the Italian Space Agency and the Netherland Space 
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Table 1. GRB Sample Chosen for the Time- Resolved Spectral Analysis 



GRB 


Redshift 


Fluence 
(xl0~ 6 cgs) 


No. of Intervals 


No. of Useful Intervals 


970111 




39.18±0.08 


10 


8 


980329 




37.53 ±0.07 


8 


6 


990123 


1.60 


205.12 ±0.03 


26 


19 


990510 


1.619 


15.80 ±0.07 


15 


8 



Note. — For each GRB we report the redshift, the fluence in the 2-2000 keV 
energy band, the number of time intervals in which we subdivided the time profile 
and those used in the time-resolved spectral analysis. 
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Table 2. Best fit Parameters of the BF in the Time Intervals for Which a Sensitive 

Spectral Analysis Was Possible 



GRB Interval Start (At) 


a 




E p 


X 2 /dof Flux 


(SOD) 






(keV) 


(x!0~ 6 cgs) 



970111 



980329 



990123 



3 


35045.0 


(3) 


^7+ - 08 
u -°' -0.06 


_o. qq+0-32 

°- yy -0.43 


2031^ 


114.2/97 


1.98 ±0.55 


4 


35048.0 


(6) 


n 90+O.O6 
u - zo -0.02 


-4.89l°:^ a 


174±| 


87.2/75 


1.79 ±0.20 


5 


35054.0 


(3) 


n go+0.05 
u - oz -0.03 


Q Q7+4.35a 
y -°' -0.63 


162+5 


89.2/73 


2.39 ±0.18 


6 


35057.0 


(5) 


-0.48lo.oj 


—4 qA+0-56a 
4 - y4 -5.06 


15113 


83.0/73 


2.44 ±0.16 


7 


35062.0 


(3) 


-0.5818:81 


_o 4fi +0.15 


89+3 


29.3/13 


1.02 ±0.14 


8 


35065.0 


(4) 


-o.58i8;8l 


_o 74 +o.2i 

a - '^-0.23 


64+2 


32.3/13 


0.58 ±0.07 


9 


35069.0 


(4) 


-o.6H8:8^ 


q 1 c+0.09 


56j* 


15.6/12 


0.71 ±0.38 


10 


35073.0 


(13) 


-0.69±g : g| 


_n opr+0.08 
z - oo -0.07 


57±| 


22.4/13 


0.43 ±0.15 



2 


13477.0 


(4) 


-0 


7/1+0.09 
'^-0.07 


9 17+O.IO 
z - 1 '-0.12 


25611 


105.7/83 


1 


50 


±0 


49 


3 


13481.0 


(2) 


-0 


c-i+0.06 
Di -0.05 


-2 46+ - 08 
z - 4D -o.io 


236l^ 9 


108.1/90 


3 


25 


±0 


66 


4 


13483.0 


(3) 


-0 


7D+0.03 
' u -0.05 


z - 4y -0.05 


235l^ 6 


105.9/97 


4 


82 


±0 


70 


5 


13486.0 


(4) 


-0 


Qn +0.05 
8U -0.02 


9 99+O.O4 
z - zz -0.04 


242l? 7 


135.5/105 


4 


25 


±0 


41 


6 


13490.0 


(4) 


-1 


9-1+0.05 
zl -0.05 


9 99+O.O7 
z - zz -0.09 


169tH 


97.7/77 





95 


±0 


15 


7 


13494.0 


(19) 


-1 


4q +0.07 
4y -0.05 


-2.79l° : ^ a 


13611 


52.8/47 





15 


±0 


04 



2 


35221 


9 


(6) 


-0 


44i8i? 


Q Q7+6.56a 
y -°' -0.63 


148+g 


39.2/35 


0.22 ±0.13 


3 


35227 


9 


(6) 


-0 


1Q+0.10 

±y -0.32 


9 ^9+0.40 a 
z - oz -7.48 


mtH 


27.8/44 


0.64 ± 0.49 


4 


35233 


9 


(2) 


-0 


cro+0.12 
b8 -0.11 


- 9 .3oi8: 7 8 a 


QSS+96 
•JOO-49 


55.1/57 


1.71 ±0.90 


5 


35235 


9 


(2) 


-0 


cc+0.03 
oo -0.03 


Q Q7+6.96a 
y -°'-0.63 


598 + | 


101.9/105 


3.80 ± 0.83 


6 


35237 


9 


(2) 


-0 


r fi +0.03 
OD -0.08 


9 f)7+0.39a 
z,.yj 1 _7 92 


8441^ 


152.8/151 


8.64 ±3.01 


7 


35239 


9 


(2) 


-0 


c:c;+0.03 
oo -0.04 


9 19+O+9 
Z - 1Z -0.71 


979l^ 7 9 


177.3/182 


13.87 ±2.60 


8 


35241 


9 


(2) 


-0 


rr/,+0.02 
°^-0.04 


9 9Q+0.28a 

— . — ■ ' y 72 


10941?^ 


164.7/183 


14.70 ± 2.91 


9 


35243 


9 


(2) 


-0 


ci+0.04 
dl -0.06 


-1 98+° :i3 
1 - yo -0.37 


5941^° 


192.0/166 


8.58 ±2.41 


10 


35245 


9 


(2) 


-0 


cc+0.02 
OD -0.03 


Q Q7+5.36a 
y .°' -0.62 


47811 


149.3/140 


5.19 ±0.68 


11 


35247 


9 


(2) 


-0 


fifi +0.07 
OD -0.04 


— 3 66 +1 ' 10a 

°. DD -6.33 


272+j° 


98.9/104 


2.12 ±0.61 


12 


35249 


9 


(2) 


-0 


r 1+0.05 
dl -0.04 


y - ou -0.63 


43111 


136.5/118 


4.35 ±0.81 


13 


35251 


9 


(2) 


-0 


44+0-05 

^-0.04 


9 90+O.I8 
Z - ZO -0.29 


5821^5 


211.4/188 


9.53 ±2.11 


14 


35253 


9 


(2) 


-0 


tro+0.04 
°°-0.03 


9 97+O.2O 
z - z '-0.35 


7-1 9+63 
' iz -69 


200.0/180 


9.84 ± 1.67 


15 


35255 


9 


(2) 


-0 


ei+0.04 
ui -0.03 


-2-441818 


6101S 


169.1/182 


7.54 ±1.26 


16 


35257 


9 


(2) 


-0 


7 i +0.04 
' 1 -0.04 


-2-77l?:^ a 


4781^0 


140.6/125 


5.13 ± 1.09 
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Table 2 — Continued 



GRB Interval Start (At) a (3 E p X 2 /dof Flux 

(SOD) (keV) (xl(r 6 cgs) 



^0.02 _o 9Q +0.67a A-\ A + U 

-0.03 °- zy -6.71 ^ 1 ^-34 

^0.03 9 nn+ - 12 9fin+ 50 

-0.03 z - uu -0.19 ^OU_ 37 

^0.03 _ 9 fi9 +0.56a qci+57 

-0.03 z - uz -7.37 ooi -55 

1-0-01 _ Q o R +6.79a / infi+ 70 

-0.03 y - ou -0.64 ^ uu -46 



990510 



17 


35259.9 


(5) 


-0. 


.91; 


18 


35264.9 


(5) 


-1 


.or 


19 


35269.9 


(5) 


-1 


.08; 


20 


35274.9 


(5) 


-1 


13; 


1 


31745.9 


(5) 


-0. 


.67; 


2 


31750.9 


(5) 


-0. 


.96; 


6 


31770.9 


(15) 


-0. 


.97; 


7 


31785.9 


(5) 


-0. 


.80; 


8 


31790.9 


(3) 


-1 


.14; 


9 


31793.9 


(4) 


-1 


.or 


10 


31797.9 


(5) 


-1 


.14; 


11 


31802.9 


(4) 


-1 


.70; 



160.5/170 


3.02 


± 


0.42 


164.4/147 


2.16 


± 


0.42 


140.8/147 


1.96 


± 


0.37 


112.9/134 


1.73 


± 


0.28 



-0.09 
-0.09 
-0.12 
-0.07 
-0.07 
-0.06 
-0.08 
-0.06 
-0.03 
-0.04 
-0.18 
-0.35 
-0.05 
-0.12 
-0.10 
-0.05 



2 


oq+0.15 
oy -0.28 


8o;? 


73.9/48 


0. 


.42 


± 


0. 


.16 


3 


cc;+0.33 
■ 00 -1.63 




53.4/53 


0. 


.54 


± 


0. 


.20 


2 


fin +0.22 
OU -0.38 




95.3/81 


1 


.29 


± 


0. 


.40 


2 


C7+O.IO 
01 -0.11 


178+H 


100.3/80 


2 


.05 


± 


0. 


,63 


9 


07+6. 31a 
°'-0.63 


151±g 


70.9/62 





.63 


± 


0. 


.13 


2. 


qc:+0.05 
OJ -0.06 


31±IS 


28.7/39 


1 


,07 


± 


1. 


,03 


2 


OR+0.11 

.86_ 12 


64+ 5 


20.0/16 





,45 


± 


0. 


.13 


9. 


07+6. 58a 
•°'-0.63 


49111 


40.9/34 





.11 




0. 


.04 



Note. — Uncertainties are single parameter errors at 90% confidence level. For each time interval 
of Figure [U we also report the Seconds Of Day (SOD) in correspondence of its start, its duration 
At in seconds, and the corresponding flux in the 2-2000 keV energy band. For the flux estimate, 
in the cases the high energy index j3 could not be constrained, it was frozen to —2.3. 



1 Value fixed to —2.3 in the computation of flux and luminosity. 
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Table 3. Correlation Analysis Results between E p and 2-2000 keV Flux 



GRB 


P 


NHP 


k 


m 


°cxt 


Xl (dof) 



970111 0.69 5.0 x 10~ 2 6.05 ± 0.34 0.68 ± 0.06 2.9 (6) 

5.91^78 0.6518:12 0.04lg;« 1.5(5) 

980329 0.49 0.32 3.25 ± 0.26 0.16 ±0.04 1-2(4) 

3.201^ o.i5i8:8? o.ooi^ 1.6(3) 

990123 0.92 2.68 x 10~ 8 5.51 ±0.28 0.53 ± 0.05 1.3(17) 

5.1618:^ o.46i8:8^ o.o 4 ^ 1.1 (ie) 

990510 0.88 3.85 xlO- 3 7.09 ± 1.01 0.81 ±0.15 1.8(6) 

5.541^ 0.5618:1 0.0918:08 1-4(5) 



Note. - A power-law relation between the two parameters is assumed: 
log-Ep — k + m log flux. In addition to the best fit parameters k and m, 
reported are the Spearman correlation coefficient p, NHP, and the best fit Xr 
with the degrees of freedom dof. The parameter estimate was performed using 
both the least-squares method and the likelihood method, in which case the 
extrinsic variance a ext between the E p ^ and flux is included in the log-likelihood 
function as a free parameter (see the text). 
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Table 4. Correlation Analysis Results between the Intrinsic (Redshift Corrected) Peak 
Energy E p i and the GRB Luminosity, Either for Individual GRBs with Known Redshift 
(990123, 990510) or When All the Data Points of These GRBs Are Taken Into Account 



GRB 


P 


NHP 


k 


m 


°"cxt 


Xl (dof) 



990123 0.92 2.7 x 10~ 8 2.64 ± 0.05 0.53 ± 0.05 1.3(17) 

2.7118;" 0.4618;" 0.04l°;& 1.1(16) 

990510 0.88 3.85 x 10" 3 2.46 ± 0.04 0.81 ±0.15 1.8(6) 

2.45l8;gg 0.5618;1 0.09l°;J° 1.4(5) 

full sample 0.94 1.57 x 10~ 13 2.51 ± 0.03 0.66 ± 0.03 1.6 (25) 

2.53±°;°1 0.63±8;gf 0.06l°;°* 1.2(24) 



Note. — A power-law relation between the two parameters is assumed: logE^j = 
k + m\ogL iso . In addition to the best fit parameters k and m, reported are the 
Spearman correlation coefficient p, NHP, and the best fit reduced xl with dof. The 
parameter estimate was performed using both the least-squares method and the 
likelihood method, in which case the extrinsic variance <j ext between the E p>i and 
flux is included in the log-likelihood function as a free parameter (see the text). 
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Fig. 1.— Light curves of GRBs 970111, 980329, 990123, and 990510, detected with Bep- 
poSAX WFC (2-28 keV) and BATSE (25-2000 keV). Also shown are the intervals in which 
the time resolved spectra were derived. 



- 22 - 





Interval number 



Interval number 




X " 10 

Z -1.5 





10 12 



Intorv.il number 



Ir.ti"'rv:il number 



Fig. 2. — Temporal evolution of the BF best fit parameters and of the null hypothesis 
probability (NHP),for each of the GRBs 970111, 980329, 990123, 990510. Top panel: low 
energy PL index a. Second panel from the top: high energy power-law index {3, in the cases in 
which it was constrained from the fit. In the other cases (not plotted) a value of /3 — —2.3 is 
assumed. Third panel from the top: peak energy E p . Fourth panel from the top: estimated 
2-2000 keV flux. Bottom panel: null hypothesis probability (NHP) for the tested function. 
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Fig. 3. — Dependence of the time resolved peak energy, obtained from the best fit of the BF 
to the joint WFC+BATSE spectra, on the 2-2000 keV flux measured in the corresponding 
interval. Also shown is the best fit power-law curve. 
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Fig. 4. — Behavior of the time-resolved low-energy photon index, obtained from the best 
fit of the BF to the joint WFC+ BATSE spectra, with the corresponding E p derived in the 
same spectral interval. 
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Fig. 5. — Intrinsic E P:i as a function of the 2-2000 keV luminosity for GRBs 990123 and 
990510. 
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Fig. 6. — Intrinsic E p ^ as a function of the 2-2000 keV isotropic luminosity L- lso , obtained by 
merging together all the available data for GRB 990123 and GRB 990510. The continuous 
line shows the best fit power-law slope (see Table HJ). 
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Fig. 7. — Compariso n of the derived E p ^ versus Li SO with the results reported by 



Ghirlanda et al.l ( 120101 ) . Left panel: our results compared with the best-fit power-law (con- 



tinuous line) and the 2a belt (dashed lines) obtained using 51 time-resolved spectra extracted 
from eight GRBs observed with Fermi. Right panel: our results compare d with the best 



fit pow er-law (continuous line) and the 2a belt (dashed lines) obtained by iGhirlanda et al. 
( 120101 ). using time-averaged spectra of 105 (pre- Fermi plus Fermi) events. 



